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Single-beam, single-target, and double-spin asymmetries for hard exclusive photon production on 
the proton ep —» e'p '7 are presented. The data were taken at Jefferson Lab using the CLAS 
detector and a longitudinally polarized 14 NH 3 target. The three asymmetries were measured in 165 
4-dimensional kinematic bins, covering the widest kinematic range ever explored simultaneously for 
beam and target-polarization observables in the valence quark region. The kinematic dependences of 
the obtained asymmetries are discussed and compared to the predictions of models of Generalized 
Parton Distributions. The measurement of three DVCS spin observables at the same kinematic 
points allows a quasi-model-independent extraction of the imaginary parts of the H and H Compton 
Form Factors, which give insight into the electric and axial charge distributions of valence quarks 
in the proton. 


I. GENERALIZED PARTON DISTRIBUTIONS 
AND DEEPLY VIRTUAL COMPTON 
SCATTERING 

It is well known that the fundamental particles that 
form hadronic matter are the quarks and the gluons, 
whose interactions are described by the Lagrangian of 
Quantum Chromo-Dynamics (QCD). However, exact 
QCD-based calculations cannot yet be performed to ex¬ 
plain the properties of hadrons in terms of their con¬ 
stituents. One has to resort to phenomenological func¬ 
tions to interpret experimental measurements in order to 
understand how QCD works at the “long distances” at 
play when partons are confined in nucleons. Typical ex¬ 
amples of such functions include form factors and parton 
distributions. Generalized Parton Distributions (GPDs), 
which unify and extend the information carried by both 
form factors and parton distributions, are nowadays the 
object of an intense effort of research, for a more com¬ 
plete understanding of the structure of the nucleon. The 
GPDs describe the correlations between the longitudinal 
momentum and the transverse position of the partons in¬ 
side the nucleon, they give access to the contribution of 
the orbital momentum of the quarks to the nucleon spin, 
and they are sensitive to the correlated qq components. 
The original articles, general reviews on GPDs and de¬ 
tails on the formalism can be found in Refs. hhz]. 

The GPDs are universal nucleon-structure functions, 
which can be accessed experimentally via the exclu¬ 
sive leptoproduction of a photon (DVCS, deeply virtual 
Compton scattering) or of a meson from the nucleon at 
high momentum transfer. More precisely, the virtuality 
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Q 2 of the photon exchanged with the nucleon by the ini¬ 
tial lepton, defined as 

Q 2 = ~(k - k') 2 , (1) 

where k and k' are the four momenta of, respectively, 
the incoming (e) and scattered electron (e / ), must be 
sufficiently large for the reaction to happen at the quark 
level. Figure [l] illustrates the leading-twist [SI process 
for DVCS, also called the “handbag mechanism”, on a 
proton target. The virtual photon interacts with one of 
the quarks of the proton, which propagates radiating a 
real photon. 


t = (p~p') 2 ( 2 ) 

is the squared four-momentum transfer between the ini¬ 
tial (p) and final proton (p 1 ). x + t; and x — £ are the lon¬ 
gitudinal momentum fractions of the quark before and 
after radiating the real photon, respectively. £ is defined 
as [H]: 


£ = Xb 


2 X b T X b ()' 


which at leading twist (—t << Q 2 ) becomes 


£ - 


x B 


2-x b ’ 


( 3 ) 

( 4 ) 


where x B is the Bjorken scaling variable 

Q 2 

XB ~ 2IWV 

M is the proton mass, and 

v = E e - E eJ . 


( 5 ) 

( 6 ) 


In the Bjorken limit, defined by high Q 2 , high v and fixed 
xb , the DVCS process can be factorized into a hard- 
scattering part ( 7 * * * § ** q —> 7 q') that can be treated pertur- 
batively, and a soft nucleon-structure part, described by 
the GPDs. At leading-order QCD and at leading twist, 
considering only quark-helicity conserving quantities and 
the quark sector, the DVCS process is described by four 
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FIG. 1. (Color online) The handbag diagram for the DVCS 
process on the proton ep —► e'p'j'. 


GPDs, H, H 1 E , E, which account for the possible combi¬ 
nations of relative orientations of nucleon spin and quark 
helicity between the initial and final state. 

The GPDs depend upon the variables x, £ and t. The 
Fourier transform, at £ = 0, of the t dependence of a GPD 
provides the spatial distribution in the transverse plane 
for partons having a longitudinal momentum fraction x 

cm. 

Model-independent sum rules link the first moment in 
x of the GPDs to the elastic form factors (FFs) [2]: 


J dxH(x,£,t) = Fi(t); J dxE(x, £, t) = F 2 {t) 

[ dxH(x,£,t) = G A (t); [ dxE(x,£,t) = G P (t),( 7) 


where F\(t) and F-^it) are the Dirac and Pauli FFs, and 
G A (t) and Gp(t) are the axial and pseudoscalar FFs. 

Among the three variables x, £ and t, which appear in 
the DVCS formalism, only two, £ and t , are experimen¬ 
tally accessible, since x appears only in the quark loop 
and is integrated over. Indeed, the DVCS amplitude is 
proportional to sums of integrals over x of the form: 



dxF(- par, £, t) 


1 

x — £ + ie 


± 


1 

x + £ — ie 


( 8 ) 


where F represents a generic GPD and the top and bot¬ 
tom signs apply, respectively, to the quark-helicity inde¬ 
pendent, or unpolarized, GPDs ( H,E ) and to the quark- 
helicity dependent, or polarized, GPDs ( H,E ). Each of 
these 4 integrals, which are called Compton Form Factors 
(CFFs), can be decomposed into their real and imaginary 
parts, as 


JteF = V 


dx 


=F 


x — i x + 


F(x,S,t) (9) 


W(£, t) = £, t) T F{-£, t)}, 


( 10 ) 


where V is Cauchy’s principal value integral. The in¬ 
formation that can be extracted from the experimental 


data at a given (£, f) point depends on the observable 
involved. When measuring observables sensitive to the 
real part of the DVCS amplitude, such as double-spin 
asymmetries, beam-charge asymmetries or unpolarized 
cross sections, the real part of the CFF, 3?eF, is accessed. 
When measuring observables sensitive to the imaginary 
part of the DVCS amplitude, such as single-spin asym¬ 
metries or cross-section differences, the imaginary part 
of the CFF, OmF, can be obtained. However, knowing 
the CFFs does not define the GPDs uniquely. A model 
input is necessary to deconvolute their x dependence. 

The DVCS process is accompanied by the Bethe- 
Heitler (BH) process, in which the final-state photon is 
radiated by the incoming or scattered electron and not 
by the nucleon itself. The BH process, which is not sen¬ 
sitive to the GPDs, is experimentally indistinguishable 
from DVCS and interferes with it. However, considering 
that the nucleon form factors are well known at small t, 
the BH process is precisely calculable. 

It is clearly a non-trivial task to measure the GPDs. 
It calls for a long-term experimental program comprising 
the measurement of different observables. Such a dedi¬ 
cated experimental program, concentrating on a proton 
target, has started worldwide in the past few years. Jef¬ 
ferson Lab (JLab) has provided the first measurement, 
in the valence region, of beam-polarized and unpolar¬ 
ized DVCS cross sections at Hall A providing a 
Q 2 -scaling test that supports the validity of the leading- 
order, leading-twist handbag mechanism starting at val¬ 
ues of Q 2 of 1-2 (GeV/c) 2 . Hall B provided pioneering 
measurements of beam [12) and target [13) spin asym¬ 
metries with the CLAS detector and afterwards ob¬ 
tained beam-spin asymmetries (BSA) over a large kine¬ 
matic range and with high statistics m- Beam-charge 
asymmetries, longitudinally and transversely polarized 
target-spin asymmetries (TSAs), as well as double-spin 
asymmetries (DSAs), have also been measured by the 
HERMES Collaboration FO . 

This paper is focused on the extraction of longitudinal 
TSAs, BSAs and DSAs for proton DVCS over a wide 
phase space using dedicated data taken at Jefferson Lab 
with the CLAS detector. 


II. DVCS ON A LONGITUDINALLY 
POLARIZED PROTON TARGET 

An analysis of DVCS observables, including the asym¬ 
metries of interest in this work, in terms of Fourier mo¬ 
ments with respect to the azimuthal angle, was carried 
out by Belitsky et al. [9], up to a twist-3 approxima¬ 
tion. These asymmetries allow one to extract separate 
components of the angular dependence of the ep —> e'j/y 
cross section, that are related to the Compton Form Fac¬ 
tors, defined in Eq. §■ The five-fold cross section for 
exclusive real photon electroproduction [9] 

da a 3 XBU 

dx B dydtd(j)dp e 16 7r 2 Q 2 \/1 + e 2 
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FIG. 2. Scheme to illustrate the definition of the angle <j>, 
formed by the leptonic and hadronic planes. 


depends on x B , t, the lepton energy fraction y = p- qi/p- 
k. with qi = k — k', e = 2and, in general, two 
azimuthal angles (p and tp e . The observables of interest 
here do not depend on ip e , the azimuthal angle of the 
scattered electron in the laboratory frame. The angle (p 
is formed by the leptonic and hadronic planes, as shown 
in Fig. [5] The charge of the electron is denoted with e and 
a represents the fine-structure constant. The amplitude 
T is the coherent sum of the DVCS and Bethe-Heitler 
amplitudes (7dvcs, 7bh): 

T 2 = |7bh| 2 + |7dvcs| 2 +21, (12) 

with the interference term defined as 


21 — 7dvcs7bh d- TdvcsTbh ■ 


(13) 


The azimuthal angular dependence of each of the three 
terms in Eq. (121 is given by [5j: 

„6 


|7bh| = 


„BH 


x bV 2 (^ + <?) 2 t'Pl(<l>)'P2 (<P) 

2 

+ ^ c® H cos rup + sf :H sin cp ], 


(14) 


|7dvcs| 2 = 


r{q 


,DVCS 
o 


2 


+ y l c n VCS cos ^ + 


y 2Q2 

s° vcs sin n(j )\} 


n—1 


(15) 


1 = 


x B yHVr{<p)V2{<p) 

s n sin n<f>\} , 


{Co +Z][ C S cos ncp + 


n—1 


(16) 


where V\ and V 2 are intermediate lepton propagators (for 
more details and definitions, see m). Among the Fourier 
coefficients cf, sf appearing in the previous expansions, 
the ones appearing in |7bh| 2 (P = BH) depend on the 
well-known electromagnetic FFs, while the ones appear¬ 
ing in |7dvcs| 2 (P = DVCS) and I (P = X) depend on 
the Compton Form Factors, the latter linearly. In the 
next sections, the sensitivity to the various CFFs of the 
DVCS spin observables presented in this paper will be 
outlined. 


II.1. Target-spin asymmetry 


The use of a longitudinally polarized (LP) target al¬ 
lows the extraction of the target-spin asymmetry A\ jl, 
where the first letter in the subscript refers to the beam 
polarization (“Unpolarized”, in this case) and the second 
to the target polarization (“Longitudinally polarized”, in 
this case), and which is given at twist-2 by: 


Ajl (<P) 


5 1,LP 


sin (f> 


pBH 

^0,unp 


+ ( C l,unp + Cqunp + ■•■) COS <p + ... 


(17) 

where the ellipses in the denominator represent smaller, 
higher-twist terms, and the unp subscript stands for “un¬ 
polarized”. The sirup coefficient s^lp, originating from 
the DVCS/BH interference term, at leading twist is pro¬ 
portional to a linear combination of the imaginary parts 
of the four CFFs (Eq. ©)> 


si.lp cx ^sm[F 1 'H + £(Fi + F 2 )(U + ^y-£) + 


(18) 


Due to the relative values of the proton form factors F\ 
and F- 2 - and given the typical JLab kinematics, the co¬ 
efficients in Eq. (18) enhance the contribution to Ajjl 
coming from SmH over the ones from other CFFs. How¬ 
ever, given that 3m 'H is expected to be about twice as 
big as 3niH, Ajjl will also be sensitive to 3 mH . In 
the kinematical range of the data presented here, higher- 
twist effects would appear in Eq. © as additional (p - 
dependent terms, the dominant of which is a sin2^> term 
in the numerator [S], 


II.2. Beam-spin asymmetry 


The expression at twist-two of the beam-spin asymme¬ 
try is 

sf,unp sirup 


Alu (<P) 


+ (C?unp + <U, unp + -) COS <p... 


(19) 


^OjUnp ' \ u l,unp ' °l,unp 


where 


4 iUnp (x 3m [F& + + F 2 )H - ^F 2 £], (20) 


Thus, the beam-spin asymmetry for DVCS on a proton 
target is particularly sensitive to the imaginary part of 
the CFF of the unpolarized GPD H. As for the TSA, 
more terms must be added if one goes beyond the leading- 
twist approximation, the larger of which is a sin 2 <p term 
in the numerator. As mentioned in Section [TJ the slope 
in t of the GPDs is related via a Fourier transform to 
the transverse position of the struck parton. Therefore, 
a measurement of 3mK and 3m FL provides, respectively, 
information on the axial and electric charge distributions 
in the nucleon as a function of x B (see Eq. 0 ) . A mea¬ 
surement of both the TSA and BSA at the same kine¬ 
matics is needed to truly distinguish between the two 
contributions. 
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II.3. Double-spin asymmetry 

The use of a polarized electron beam along with a lon¬ 
gitudinally polarized target allows also the determina¬ 
tion of the double-spin asymmetry All- Unlike Aul, the 
Bethe-Heitler process alone can generate a double-spin 
asymmetry. At twist-2 level, this observable takes the 
form: 


A L l(0)■ 


BH , / BH , k 

1 0,unp ' V^yunp ' °l,unp 




with 


cJ,lp> 'u.lp of SJe/Fj’H + £(Fi + F^iji H —— £) + 




( 22 ) 


where terms depending on powers of were neglected. 
In this expression, the interference terms, related to the 
GPDs, are expected to be smaller than the known BH 
terms 0. Moreover, both the constant and the cos in¬ 
dependent terms contain contributions from both BH and 
the DVCS/BH interference. Nonetheless, given the fast 
variation of BH depending on the kinematics, it is im¬ 
portant to sample All over a wide phase space to find 
possible regions of sensitivity to 3 ?eH and 3 fteV., which 
dominate in Eq. (22). 


III. THE EXPERIMENT 

The data were taken in Hall B at Jefferson Lab from 
February to September 2009, for a total of 129 days. A 
continuous polarized electron beam was delivered by the 
CEBAF accelerator onto a solid ammonia target, polar¬ 
ized along the beam direction. Frozen beads of param- 
agnetically doped 14 NH 3 , kept at temperatures of about 
IK and in a 5T magnetic field, were dynamically po¬ 
larized by microwave irradiation HU. The target was 
1.45 cm long and 1.5 cm in diameter. The target system 
included a 4 He evaporation refrigerator and a supercon¬ 
ducting split-coil magnet. The magnet, which provided 
a uniform polarizing field (A B/B = 10“ 4 ) and at the 
same time focused the low-energy Mpller electrons to¬ 
ward the beam line, was inserted in the center of the 
CLAS detector m ■ CLAS was a spectrometer, based 
on a toroidal magnetic field, providing a wide accep¬ 
tance. The CLAS magnetic field was generated by six 
superconducting coils arranged around the beamline to 
produce a field pointing primarily in the azimuthal di¬ 
rection. The particle detection system consisted of drift 
chambers (DC) to determine the trajectories and the mo¬ 
menta of charged particles curved by the magnetic field, 
gas Cerenkov counters (CC) for electron identification, 
scintillation counters for measuring time-of-flight (TOF) 
and electromagnetic calorimeters (EC) to detect shower¬ 
ing particles (electrons and photons) and neutrons. The 


segments were individually instrumented to form six in¬ 
dependent magnetic spectrometers with a common tar¬ 
get, trigger, and data-acquisition system. An additional 
detector, the Inner Calorimeter, constructed for a pre¬ 
vious DVCS-dedicated experiment m to complete the 
photon acceptance at low polar angles (from 4° to 15°), 
was placed at the center of CLAS. Figure [3] shows the 
whole setup, including the polarized target, CLAS and 
the IC. A totally absorbing Faraday cup (FC), down¬ 
stream of CLAS, was used to determine the integrated 
beam charge passing through the target. 

The data presented here come from the first two of the 
three parts in which the experiment was divided, and for 
which a 14 NH 3 target was used. The main differences 
between these two parts (referred to as A and B) were 
the target position with respect to the center of CLAS 
(z = —57.95 cm for part A and 2 = —67.97 cm for part 
B) and the electron-beam energy (F = 5.886±0.005 GeV 
for part A and E = 5.952 ± 0.005 GeV for part B). The 
beam energy values were obtained from elastic-scattering 
analysis on these data, using the 14 NH 3 target. 

The beam was rastered over the target in a spiral mo¬ 
tion in order to assure a homogeneous depolarization 
over the whole volume of the target. The beam po¬ 
larization was frequently monitored in Mpller runs, via 
the measurement of the asymmetry of elastic electron- 
electron scattering. The target polarization was con¬ 
tinuously monitored by a Nuclear Magnetic Resonance 
(NMR) system. Runs on a carbon target were taken 
periodically throughout the duration of the experiment 
for unpolarized-background studies. The selection of the 
data to be analyzed was done by monitoring the stabil¬ 
ity over time of Faraday-cup normalized yields per sector. 
For the results presented here, runs taken with the 14 NH 3 
and carbon targets from parts A and B were used. The 
two parts were analyzed separately and the final results 
were found by combining the results from the two parts, 
as will be described in Section IV.91 


IV. DEFINITIONS OF THE ASYMMETRIES 


This paper reports on the extraction of three kinds of 
asymmetries, the experimental definitions of which are 
given here. In all of the formulae below, the first sign in 
the superscript on the number of normalized DVCS/BH 
events N is the beam helicity (6) and the second sign 
is the target polarization (t). N is obtained from epy 
events (V ep7 ), normalized by the corresponding Faraday- 
cup charge ( FC bt ) after subtraction of the 7r° background 
as follows: 


N bt 


(1 


N bt 

ep'Y 

FC bt ’ 


(23) 


where B^o is the relative 7r° contamination, outlined 
in Section IV. 71 

The beam-spin asymmetry is calculated as: 
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FIG. 3. Drawing of the experimental setup, including the CLAS detector with its components (DC, EC, CC, TOF), the Inner 
Calorimeter and the polarized target. 


A P t ~(N++-N-+) + P t + (N+- -N-) 

LU P b (P t ~(N++ +N~+) + P t + (N+~ + N )) ’ 

_ (24) 

where P b is the polarization of the beam and P f +< ' ^ are 
the values of the polarization of the target for its two 
polarities. 

The target-spin asymmetry is computed as: 

^4-ul = ^ul + c a ut j (25) 


The double (beam-target) spin asymmetry is obtained 
as: 


All = A 


lab 

LL 


CAr, 


(27) 


where 


where 


a lab _ 

^UL — 


N++ + N~+ - N+- - N~ 


D f (P t ~{N++ + N~+) + P t + {N+~ + N )) 


(28) 

Df is the dilution factor to account for the contribution 
of the unpolarized background (Section V.5) and ca vt 
represents a correction applied to define the target-spin 
asymmetry with respect to the virtual-photon direction 
(Section V.lOl. 


a lab _ 

^LL — 


N++ + N— - N+~ - N~+ 


P b ■ D f (P t ~{N++ + N~+) + P t + {N +- + TV—)) 

(28) 

and ca lt is the analo g of ca ut for the double-spin asym¬ 
metry (Section V.lOl. 

In the following, the steps leading to the extraction from 
the data of all the terms composing these asymmetries 
will be presented. 
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V. DATA ANALYSIS 


V.l. Particle identification 


The final state was selected requiring the detection of 
exactly one electron and one proton, and at least one 
photon. 

The electrons were selected among all the negative 
tracks with momenta above 0.8 GeV/c, but requiring 
that their energy deposited in the inner layers of the EC 
■ IS] be greater than 0.06 GeV, in order to reject nega¬ 
tive pions, and that their hits in the CC and in the SC 
be matched in time. Fiducial cuts were also applied to 
eliminate the events at the edges of the EC (where the 
energy of a particle cannot be correctly reconstructed be¬ 
cause a large part of the induced electromagnetic shower 
is lost), and to remove the “shadow” of the IC, which 
limits the CLAS acceptance for charged particles. The 
effect of these cuts on the distribution of the total en¬ 
ergy deposited in the EC divided by the momentum is 
shown as a function of the momentum in Fig. HI The 
comparison of the top and middle plots of Fig. [5] shows 
the effect of these same cuts on the number of CC photo¬ 
electrons. Most of the events in the single-photoelectron 
peak, which come from either pions or electronic noise 
in the PMTs, are removed by our particle identification 
(PID) cuts for electrons. The rest of the background is 
eliminated by the e xclusi vity cuts applied afterwards as 
explained in Section V.3 (Fig. [5| bottom). 

The protons were selected from the correlation between 
their velocity, deduced from the time-of-flight measure¬ 
ment using the SC, and the proton momentum as mea¬ 
sured by the DCs (Fig. [6]). Fiducial cuts on 9 and were 
also applied in order to remove the shadow of the IC. 

For the photons, two different sets of cuts were 
adopted, depending on whether the photon was detected 
in the IC or in the EC. For the IC case, a low-energy 
threshold of 0.25 GeV and a cut on 9 versus E 1 to remove 
the background coming from Mpllcr electrons (Fig. [7]) 
were applied, as well as fiducial cuts on x and y , to re¬ 
move the outer and inner edges of the detector, where 
clusters could not be fully reconstructed. For the EC 
case, all neutrals passing a low-energy threshold of 0.25 
GeV and having /3 > 0.92 (Fig. [8]) to select the in-time 
photons were retained. Fiducial cuts as for the electrons 
were also adopted to remove the edges of the detector 
and IC-frame cuts were applied to remove the shadow of 
the IC over the EC. 


V.2. Monte-Carlo simulations 



p (GeV/c) 


FIG. 4. (Color online) Total energy deposited in the EC (in- 
ner+outer layers), E, divided by the particle momentum p 
as a function of p for all the negative tracks. Top: nega¬ 
tive charged particles, before cuts. Bottom: after minimum- 
momentum, ECi„, fiducial and timing cuts. 


• optimization of the DVCS and epir 0 exclusivity cuts 
(Sections V.3 and V.7 respectively); 


evaluation of the ep7r° contamination in the epj 


event sample (Section V.71. 


The two sets of gene rated e vents, DVCS/BH and epn° 
(Section V.2.1 and V.2.2), were fed to the CLAS 


GEANT3-based simulation package (GSIM) to simulate 
the response of the CLAS detector. The output of GSIM 
was then fed to a post-processing code (GPP) that simu¬ 
lates dead or inefficient DC wires and smears the DC and 
TOF resolutions to more closely agree with the experi¬ 
ment. The output of GPP was finally fed to the CLAS re¬ 
construction package, and reduced ntuples and root files 
were obtained from the reconstructed files in the same 
way as was done in the data processing. 


Monte-Carlo simulations were used in this analysis for 
the three following tasks: 

• determination of angular and momentum correc¬ 
tions for electrons and protons to compensate for 
multiple scattering and energy losses in the target 
and detector materials, as well as for imperfections 
in the trajectory reconstruction; 


V.2.1. DVCS/BH simulation 

A DVCS/BH event-generator code, which produced 
epq events according to the formalism of Belitsky et al. 
0, was used for the Monte-Carlo simulation. Figure 
[9] shows a comparison of the distributions of the rele¬ 
vant kinematic variables for the data (black lines) and 
the Monte-Carlo simulation (shaded areas). Here, PID 



















FIG. 5. Number of CC photoelectrons (times 10) for all neg¬ 
ative tracks (top), after applying all PID and fiducial cuts for 
electrons (middle), and after epy exclusivity cuts. 



p (GeV/c) 


FIG. 6. (Color online) (3 as a function of p for positively 
charged particles. The black lines represent the cut applied 
to select protons. 



FIG. 7. (Color online) Polar angle 9 as a function of the 
reconstructed energy for IC hits, showing the cut on the min¬ 
imum energy at 0.25 GeV, as well as the “triangular” cut to 
remove the low-energy/low-6 1 accidental background, applied 
to select photons. 



Energy (GeV) 


FIG. 8. (Color online) Distribution of (3 vs energy for neutrals 
measured by time-of-flight with EC. The events for which 
/3 > 0.92 (black line) were retained as photon candidates. 


and ep7 exclusivity cuts, which will be described in Sec¬ 
tion |V.3[ were included for both the data and the Monte 
Carlo. The agreement between data and simulation is 
quite good, especially given the purposes of the simula¬ 
tions in these analyses: they are not used for absolute 
acceptance corrections, but only to help in the deter¬ 
mination of cut widths and for background subtraction. 
Slight differences between data and MC, especially visi¬ 
ble in the high-t and central-^ region, are coming from 
events in the EC topology. This is probably due to the 
fact that the data, unlike the MC, are not only pure 
epy events, but are contaminated by exclusive n 0 events. 
The fact that the ep7r° contamination is larger in the EC 
topology, as will be reiterated in Section [VT3j can explain 
the data/MC discrepancies. 

V.2.2. Exclusive tt° simulation 

Exclusive epn 0 events were generated using a code for 
meson electroproduction that included a parametrization 
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x b <|> (deg) 

FIG. 9. DVCS channel. Comparison of data (black lines) and 
Monte-Carlo simulation (shaded areas). From the top left, 
Q 2 , —t, xb and <j> are plotted. The histograms are normalized 
to each other via the ratio of their integrals. 


of the ep7r° differential cross sections that have recently 
been measured by CLAS 119j . Figure [l0| shows a compar¬ 
ison of the distributions of the relevant kinematic vari¬ 
ables for the data (black) and the Monte-Carlo simula¬ 
tion (shaded areas). Here, PID and ep7r°-exclusivity cuts, 
which will be described in Section |V.7[ were included for 
both the data and the Monte Carlo. The agreement be¬ 
tween data and simulation is quite good. 

V.3. DVCS channel selection 

After selecting events with exactly one electron and 
one proton, and at least one photon, and having ap¬ 
plied the momentum corrections, further cuts need to be 
applied to ensure the exclusivity of the epy final state. 
Two kinds of backgrounds need to be minimized: the 
nuclear background coming from scattering off the ni¬ 
trogen of the 14 NH 3 target, and the background coming 
from other channels containing an electron, a proton and 
at least one photon in the final state. Having measured 
the four-momenta of the three final-state particles, one 
can construct several observables, hereafter referred to 
as “exclusivity variables”, on which cuts can be applied 
to select the epj channel. Here, the following quantities 
were studied: 

• the squared missing mass of the ep system, 
MM 2 (ep); 

• the angle 0 7 x between the measured photon and 
the calculated photon, using the detected electron 



x b <i> (deg) 

FIG. 10. Exclusive n° channel. Comparison of data (black 
lines) and Monte-Carlo simulation (shaded areas). Starting 
from the top left, Q 2 , — t, xb and (j> are plotted. The his¬ 
tograms are normalized to each other via the ratio of their 
integrals. 


and proton together with ep —> e'p'X kinematics; 

• A <j>, the difference between two ways to compute 
the angle </> between the leptonic and the hadronic 
plane. The two ways concern the definition of the 
normal vector to the hadronic plane: one is via the 
cross product of the momentum vectors of the pro¬ 
ton and the real photon, and the other one is via the 
cross product of the momentum vectors of the pro¬ 
ton and the virtual photon. For the epy final state, 
A (j> should be distributed as a gaussian centered at 
zero, with width determined by the experimental 
resolution; 

• Pperp , the transverse component of the missing 
momentum of the reaction ep —> e'p'^X, given 
by Pperp = \/Px{ X) 2 +p y (X) 2 , in the laboratory 
frame. 

The definition of the exclusivity cuts is quite a delicate 
step in the DVCS analysis. It is important that the cuts 
are determined in a consistent way for the data and for 
the Monte-Carlo simulation, because the latter will be 
used to evaluate the n° background contamination. In 
this analysis it was chosen to define each exclusivity cut 
by fitting the corresponding variable with a gaussian and 
cutting at ±3<r around the fitted mean. This procedure 
was done separately for the data and for the simulation. 
This way, the same fraction of events was kept, for both 
data and simulation. The exclusivity variables to be fit¬ 
ted were plotted after applying preliminary cuts that in¬ 
cluded: 
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• Kinematic cuts to be above the region of the nu¬ 
cleon resonances and to approach the regime of 
applicability of the leading-twist GPD formalism: 
Q 2 > 1 (GeV/ c) 2 , -t < Q 2 , an d W > 2 GeV/c 2 ; 
(where W = \JM 2 + 2 Mu — Q 2 ) 


• Ej > 1 GeV, since the real photons of interest are 
expected to have high energy; 


Bin 

xb bin 

9 e bin 

(x B ) 

(O') ((GeW) 

1 

0.1 <x B < 0.2 

15“ 

< 9 e < 48° 

0.179 

1.52 

2 

0.2 < xb < 0.3 

15° 

<9 e < 34° 

0.255 

1.97 

3 

0.2 <x B < 0.3 

34° 

<9 e < 48° 

0.255 

2.41 

4 

0.3 < x B < 0.4 

15° 

< 9 e < 45° 

0.345 

2.60 

5 

xb > 0.4 

15° 

<9 e < 45° 

0.453 

3.31 


TABLE I. Definition of the bins in xb and 9 e (Q 2 ), and av¬ 
erage kinematics for xb and Q 2 ((GeV/c) 2 ) for each bin. 


• 3cr cut around the mean of MM 2 (ep) to eliminate 
from the experimental data the background from 
channels other than epy or ep7r° (visible in the top 
left plot of Figs. [TT] and [T2| where peaks from p and 
ut/p are evident). 


In order to eliminate broadening on the widths of the 
peaks due to events from electron scattering on the ni¬ 
trogen, the fits to the exclusivity variables were done on 
the spectra obtained after subtracting carbon data from 
the 14 NH 3 data. The two datasets were normalized to 
each other via the ratio of their Faraday-cup counts mul¬ 
tiplied by a constant that accounts for different densities 
of materials for the two target types (see Section V.5). 


The method to define the cuts described above was 
adopted for the topology where the photon was detected 
in the IC, since the comparison with Monte Carlo showed 
that these data are strongly dominated by the DVCS/BH 
channel. Figure fTT| shows, for the IC topology, the effects 
of the exclusivity cuts, which appear successful both in 
extracting quite a clean epy final state (shaded areas) 
and in minimizing the background originating from the 
nitrogen part of the target (black areas). 


A different strategy was found to be necessary for the 
EC case, which displayed, before cuts, a larger contri¬ 
bution from epn° events. The peaks in the exclusivity 
variables for the data in this topology are very broad, 
when visible, and not necessarily produced by DVCS/BH 
candidates. In the first plot of Fig. [l2j for example, the 
distribution of the squared missing mass of the ep system 
is shown for the EC case. As is clearly visible, the peak 
of the distribution is not centered at zero but around the 
squared 7r° mass, indicating a significant contamination 
from the exclusive 7r° events that will be subtracted later 
through the procedure described in Section |VTj 

Thus it was decided, for the EC topology, not to fit 
the distributions of the exclusivity variables to extract 
cut means and widths. Instead it was chosen to fit only 
the peaks of the DVCS/BH Monte-Carlo simulations. To 
correct for the discrepancies in resolutions between data 
and simulation, the widths of the various exclusivity vari¬ 
ables obtained from the fits were then multiplied by ap¬ 
propriate scaling factors. These factors were obtained 
from the comparisons of data/MC for the epit 0 channel 
in the EC-EC topology. The cuts and their effects are 
shown in Fig. [I2| 


V.4. Four-dimensional binning and central 
kinematics 


The DVCS reaction can be described by four indepen¬ 
dent kinematic variables. The typical variables used to 
interpret the results in terms of Generalized Parton Dis¬ 
tributions are Q 2 , Xb, —t and </. In accordance with the 
choice made in previous DVCS analyses HE the binning 
of the data in the Q 2 -x b plane was done making 5 slices 
in the polar angle 9 e of the electron and in xb ■ The limits 
of the slices are given in Table [TJ as well as the bin aver¬ 
ages for Q 2 and xb , defined as the weighted average over 
the distribution of events in each bin. The size of the 
bins was optimi zed to have comparable statistics. The 
top plot of Fig. 13 shows the chosen grid in the Q 2 -Xb 
plane. Ten equally spaced bins in <j> and 4 bins in —t were 
adopted. The bin limits and data-averaged bin centers 
are summarized in Table 0 The bottom plot of Fig. [I3| 
shows the binning in the t-xs plane. 


Bin 

—t range (GeV/c)^ 

(- 1 ) (GeV/c) 2 

1 

0.08 < -t < 0.18 

0.137 

2 

0.18 <-t< 0.3 

0.234 

3 

0.3 < -t < 0.7 

0.467 

4 

0.7 < -t < 2.0 

1.175 


TABLE II. Definition of the bins in —t and average kinematics 
for each bin. 


The central kinematics in this analysis were defined as 
the average value from the data of each of the 4 kinematic 
variables for each bin. In fact, at first order the uncertain¬ 
ties on an asymmetry induced by taking bins of finite size 
are minimized when the central kinematics are chosen to 
be the weighted average over the distribution of events 
in that bin. The procedure to compute second-order bin¬ 
centering corrections is reported in Section V.ll The 
grid of bins was applied to both parts A and B of the 
experiment. In order to establish whether the asymme¬ 
tries obtained from the two sets of data could eventually 
be combined, the central values of the various bins into 
which the available phase space was divided were com¬ 
puted, and compared for the two parts. For each bin in 
Q 2 (0 e )-XB and —t, the central kinematics for parts A 
and B were found compatible well within their standard 
deviations. 
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FIG. 11. IC topology. Effects of the exclusivity cuts. Top left: squared missing mass of the ep system; top right: the angle 
9 y x between the measured photon and the calculated photon from ep —> e'p'X- bottom left: difference between two ways to 
compute the (f> angle; bottom right: missing transverse momentum p per p calculated from ep —> e'p'^X. The dot-dashed and 
solid lines show the events before exclusivity cuts for, respectively, 14 NH 3 and 12 C data, while the gray and black shaded areas 
represent the events after all exclusivity cuts except for the one on the plotted variable for, respectively, 14 NH 3 and 12 C data. 
The lines and arrows show the limits of the selection cuts. The plots for 14 NH 3 and 12 C data are normalized to each other via 
their relative luminosities. 


V.5. Dilution factors 


For both the DVCS and epn° final states, dilution fac¬ 
tors to correct the experimental yields for the contribu¬ 
tion from the scattering on the unpolarized nitrogen of 
14 NH 3 were determined using data taken on 14 NH 3 and 
on 12 C targets. The dilution factor is defined as 


D f = 1 - c ’ 


Nl2Q 
Nn nh 3 


(29) 


Here, Nhq is the number of events, normalized by the 
corresponding Faraday-cup counts, taken on carbon and 
surviving all of the DVCS (or epn°) selection cuts, while 
7 V i4 nh 3 is the number of events, normalized by the cor¬ 
responding Faraday-cup counts, taken on 14 NH 3 passing 
the DVCS (or epn°) selection cuts. The factor c accounts 
for the different luminosities of the two sets of data, which 
are in turn related to the ratio of the areal densities of 
the materials present at the target level for the two kinds 
of runs ( 14 NH 3 in the numerator, 12 C in the denomina¬ 
tor). As the statistics acquired on carbon during the 
experiment was much smaller than for the 14 NH 3 data, 
it was not possible to perform the dilution factor anal¬ 
yses for each 4-dimensional bin. The dependence of the 
dilution factor on each of the 4 kinematic variables was 
checked by integrating over the other three. The results 
for part B, as an exam ple, are shown as functions of xb, 
Q 2 , —t, and tf> in Fig. 14 Both dilution factors show 


an approximately flat dependence in each of the 4 kine¬ 
matic variables. There may be a small Xb dependence, 
although the dilution factor is not inconsistent with a 
constant behavior. The fit results for each variable are 
consistent with each other within error bars. Therefore, 
a constant value of Df for all the kinematic bins, for both 
the epy and epir 0 analysis, was assumed. The following 
values were adopted: 

• part A: D/(epy) = 0.912 ± 0.009; 
part B: Df{ep r f) = 0.928 ± 0.006; 

• part A: D/(ep7r°) = 0.921 ± 0.016; 
part B Df(epn°) = 0.896 ± 0.010. 


V.6. Beam and target polarization 


The product of the polarizations of beam and target, 
PbPt , was extracted for part A and part B separately. 
Exclusive elastic events were used for this analysis, since 
the elastic asymmetry is well known [20] and PbPt is the 
proportionality factor between the experimental asym¬ 
metry ( A exp ) and the expected one (A): 


A = 


PbPt ’ 


(30) 
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FIG. 12. EC topology. See caption of Fig. [TT| 



FIG. 13. (Color online) Grid showing the binning in the Q 2 - 
xb space (top) and in the —t-XB space (bottom). 
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FIG. 14. Dilution factor as a function of (from left to right) 
xb , Q 2 , —t, <j>, for the DVCS analysis (top) and for the epiv 0 
analysis (bottom). Part B. 


where 


4 N + — N~ 

exp ~ D f (N++N ~) 


N+(~) refers to the number of elastic events for posi¬ 
tive (negative) beam helicity for a given target polar¬ 
ization, normalized to the corresponding Faraday-cup 
counts. Df is the dilution factor (~ 98%) determined 
in the same way as for the DVCS and ep —> e'p'n 0 
events (Section |V.5| ). The expected asymmetry A is com¬ 
puted according to [20], and for the elastic form factors 
a parametrization, obtained by fitting the results from 
a JLab polarization-transfer experiment [21] . was used. 
Inverting Eq. (30) PbPt was calculated for each run part 
as the average over 7 different Q 2 bins. The beam polar¬ 
ization alone (Pj,) was obtained from the analysis of the 
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Run part 

P b Pt + 

P b P t ~ 

P b 

Pt + 

Pf 

Part A 
Part B 

0.648 ±0.018 
0.674 ±0.011 

0.625 ±0.016 
0.625 ±0.010 

0.87 ±0.04 
0.837 ±0.017 

0.75 ±0.04 
0.81 ±0.03 

0.72 ± 0.04 
0.75 ±0.03 


TABLE III. Obtained values of PbPt, Pb, and P t , for each part of the experiment, and for both positive and negative target 
polarizations. 


various Mpller runs that were taken during the experi¬ 
ment. Event-weighted averages of Pb were computed for 
the two parts of the experiment. The target polarization 
values were thus deduced by dividing the measured PbPt 
by Pb- The results for the product of the polarizations 
of beam and target for each polarization sign, and of the 
beam and target polarizations alone for the two parts of 
the experiment, are summarized with their uncertainties 

in Table Hill 


V.7. n° background subtraction 

The procedure adopted to compute the 7 r° contamina¬ 
tion to the epy event sample is described here. All the fol¬ 
lowing was done separately for each of the 4 -dim ensional 
(Q 2 , Xb , —t, </>) bins described in Section V.4| and for 
each beam and target polarization state. To keep the 
notation simpler, in the following formulae we omit the 
dependence on (Q 2 , Xb , — t, 4>) of each quantity. After 
applying the epy event selection cuts and the DVCS ex¬ 
clusivity cuts, one is left with N epl events that are not 
only DVCS or Bethe-Heitler candidate events, but can 
also contain some epn 0 events in which one of the two 
7 r°-decay photons has escaped detection. The number of 
DVCS/BH events can be obtained as 

-^dvcs/bh = N epl (\ - B n o), (32) 

where the contamination fraction is given by 
D f (epTT°)N ep ^o (l7) 


B~ o = 


Df{em)N e . 


(33) 


To calculate N epn o( l7 ), we use both the MC simulation 
and exclusive epn 0 events {N® A ^ A ) from the real data: 


N epn o ( i 7) = N" A1 0 A * i?Acc(l7/27)- 


(34) 


The ratio of acceptances -Ra<x( 17/27) is obtained from 
simulated epn 0 events, and is defined as: 


n mg 

D /o..\_ ep7T°(I7) 

uAcc(l7/27) — t/MC ■ 

ep-7r°(27) 


(35) 


where is the number of reconstructed events ob¬ 

tained by applying the DVCS selection cuts to the epn° 
simulation, and N^^ 0 r 2 -y) the number of events ob¬ 
tained by applying the epir 0 selection cuts to the same 
epn 0 simulation and for the same number of generated 
events. 


^Vep 7 r° A was extracted by selecting events with one elec¬ 
tron, one proton and at least two photons using the same 
PID cuts as for the DVCS channel. The chosen photon 
pair would then be the one whose invariant mass was the 
closest to the nominal 7 r° mass. Then, the same method 
as the o ne a dopted to define the DVCS selection cuts 
(Section V.3| was adopted to obtain the exclusivity of 
the epn 0 final state. Three “exclusivity variables” were 
used: 


• the missing mass squared of the ev system, 
MM 2 (ep); 

• the angle between the measured and calculated 7r° 
direction, 9^o X ', 

• the two-photon invariant mass, 1 X 1 ( 77 ). 

The carbon data were subtracted from the 14 NH 3 data 
to remove possible smearing effects on the fitted peaks 
coming from background events due to electron scatter¬ 
ing on 14 N. When possible, data and Monte-Carlo distri¬ 
butions were fitted separately, at the “preliminary cuts” 
stage (the preliminary cuts being the same as for the 
DVCS selection, excluding the E 7 > 1 GeV cut). The 
obtained widths were compared to the widths of the dis¬ 
tributions after all cuts were applied to check for possible 
correlation effects. The three different photon-detection 
topologies were treated differently, as in the case of the 
DVCS selection. Figures [15] and [16] show the results of 
the epn 0 exclusivity cuts for the data when both photons 
were detected in the IC and in the EC, respectively. 

The acceptance ratio Race varies for each kinematic 
bin, and is, on average, around 30%. In absolute value, 
the background/signal ratios range from 1 % (at low t, 
low-mid xb, high Q 2 ) to about 50% (in the highest t, 
highest xb, central <f> range). However, what will im¬ 
pact the final asymmetries is not the size of the contam¬ 
ination itself, but the point-by-point difference of con¬ 
tamination for positive and negative target (or beam, or 
beam-target) polarization. The background subtraction 
was done for each individual 4-dimensional bin. The av¬ 
erage impact of the background subtraction on each of 
the 3 asymmetries, relative to their average value at 90°, 
is ~ 8 % for the BSA, ~ 11% for the TSA, and ~ 6 % for 
the DSA. 


V.8. Nuclear background effects on BSA 


The beam-spin asymmetry extracted using the 14 NH 3 
data contains a contribution from the nitrogen nuclei. 
The dilution factor analysis (Section V.5| shows that, 
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FIG. 15. Exclusive tt° analysis, IC-IC topology. Effects of 
the epn° exclusivity cuts on MM 2 (ep) (top), 6 ^ 0 x (middle), 
and two-photons invariant mass (bottom). The dot-dashed 
and solid lines show the events before exclusivity cuts for, 
respectively, 14 NHs and 12 C data, while the gray and black 
shaded plots represent the events after all exclusivity cuts but 
the one on the plotted variable for, respectively, 14 NH 3 and 
12 C data. The lines and arrows show the limits of the selection 
cuts. 


thanks to the tight exclusivity cuts adopted, less than 
10% of the final DVCS/BH event candidates come from 
the nitrogen. However, given the low statistics collected 
in the carbon runs, it is not possible to directly subtract 
for each kinematic bin and helicity state the nitrogen 
contribution, calculated on the basis of the 12 C measure¬ 
ment, from the 14 NH 3 events. The raw BSA, without 7r° 
subtraction, was extracted from the carbon-target data, 
applying all of the DVCS/BH selection cuts and integrat¬ 
ing, for statistics reasons, over all kinematic variables but 
cj). The resulting asymmetry was compared to the corre¬ 
sponding observable obtained on the 14 NH 3 target. A 
Student’s t-test was performed, and it was found that 
the two DVCS/BH BSAs for a free proton and quasi-free 
proton in carbon are compatible at the 3 a level. This 
result points to the fact that given the limits imposed by 


FIG. 16. Exclusive 7 r° analysis, EC-EC topology. See caption 
of Fig. [15] 


our statistics, our experiment is not sensitive to possible 
medium-modification effects for DVCS, going from pro¬ 
ton to carbon target. Moreover, our results for the dilu¬ 
tion factor show a contamination of events from nitrogen 
in the epy sample of less than 10%, and this contamina¬ 
tion comes from bound protons at quasi free kinematics 
due to the tight exclusive cuts. All this supports the 
assumption that the contribution to the final beam-spin 
asymmetries of the incoherent DVCS on nitrogen does 
not impact sizeably our result, within the limits of our 
statistical accuracy. 


V.9. Combination of asymmetries from part A and 
part B 


Having verified that the bin centers for parts A and 
B are equal within uncertainties (Section V.4), the final 
asymmetries, after epn 0 background subtraction, have 
been combined via uncertainty-weighted averages for 
each four-dimensional bin in (Q 2 , Xb, —t, <j>). In order 
to verify that the two sets of asymmetries are statisti- 
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cally compatible with each other, f-tests were perfomed 
for the three asymmetries. The tests proved that the 
fluctuations between part A and B are purely statistical, 
and the two datasets can therefore be combined. The 


combined asymmetries are shown in Section VII 


The combined beam energy, 5.932 GeV, was computed 
as the weighted average over all epy events of the values 
of the energies for the two parts, which were separately 
deduced from elastic measurements (Section |III|). 


V.10. Transverse corrections 


In this experiment, the target polarization was par¬ 
allel to the electron beam direction. However, for the 
theoretical interpretations of the asymmetries, the lon¬ 
gitudinal polarization with respect to the virtual photon 
direction is usually adopted. In order to be consistent 
in the comparisons with the theoretical models, a model- 
dependent correction was computed to finally obtain the 
TSA (and DSA) with respect to the virtual photon di¬ 
rection. The purely longitudinal asymmetry referred to 
the virtual photon, Aul, is linked to the ^-component of 
the transverse one, Aut(O), and to what is measured in 
the lab, Ay}), by the relation [22] : 

A lab 

Aul = —+ tan0*A UT (O), (36) 

cos 9* 

where 9* is the angle formed by the virtual photon and 
the beam direction. An analogous relation holds for the 
double-spin asymmetry, 

Alab 

All = —+ tan0*A L T(O). (37) 

costr 


The idea adopted here is to compute a set of model- 
dependent bin-by-bin corrections using the predictions 
of the GPD model VGG @2] (described in Section VI) 


for Ajjt{ 0), and the average value of 9* over the selected 
DVCS events for each bin. The VGG model was run for 
all kinematic bins to obtain Aj/t(0), with three different 
sets of options for systematic studies of the model depen¬ 
dence of the correction. The size of the transverse-target 
spin asymmetry varies with —t, from a few percent at low 
—i to up to 30% for the highest —t bins. Also the dif¬ 
ference between the results obtained for the three sets of 
VGG options is more sizeable at high —t, where changes 
of sign for the <f> distribution are observed. As far as 
the double-transverse spin asymmetry is concerned, it is 
higher at low —t, mostly positive, and never greater than 
30%. The angle 9* between the virtual photon and the 
beam direction is, on average, around 7°. Thanks to the 
combination of small values for 9* and relatively small 
values for Aut(0), the corrections, which are defined as 


c a ut — Aul — Ay]) 

ca lt = ^ll - All 


(38) 

(39) 


three versions of “corrected” TSAs (or DSAs) and the 
measured one was adopted as a systematic uncertainty. 
In general, the correction brings a slight increase to both 
kinds of asymmetries, and it is always smaller than the 
statistical uncertainties. The same holds for the asso¬ 
ciated systematic uncertainty, the averages of which for 
both TSA and DSA are reported in Table [TV] The trans¬ 
verse correction was applied to each 4-dimensional bin. 
The average impact of this correction on each of the 2 
asymmetries, relative to their average value at 90°, is 
~ 4% for the TSA and ~ 2% for the DSA. The values of 
the corrections for each 4-dimensional bin are reported 
in Table [V] along with the values of the corrected asym¬ 
metries. 


V.ll. Bin-centering corrections 


Second-order bin-centering corrections were applied to 
the single-spin asymmetries. As stated in Section |V.4[ 
the asymmetries are correctly defined at first order if one 
takes the average value of each kinematic variable over 
all the events of each bin. To apply second-order bin¬ 
centering corrections, one needs to make assumptions on 
the kinematic dependences of the asymmetry. Based on 
the information from the present data, only the effect of 
the —t dependence was studied, since the asymmetries 
do n ot ex hibit strong variations in Q 2 and xb (see Sec¬ 
tion VIII. The bin-centering corrections were evaluated 


for each 4-dimensional bin as 


cbcc = 


a(t) sin^i 


a((t))sin(0) 


1 + f3(t) cosfi 1 + f3((t)) cos(0) ’ 


(40) 


where the first term is computed event-by-event while the 
second term is evaluated at the central values ( t ) and (</>) 
of the bin. The a(t) and (3(t ) functions were determined 
by fitting the data. 

For the BSA and the TSA, the corrections are on av¬ 
erage at the level of 4%, relative to the average value of 
the asymmetries at 90°. Since the extrema in the DSA 
are much less pronounced, the corrections are expected 
to be even smaller and they are therefore neglected here. 


V.12. Systematic uncertainties 

Systematic checks were performed to evaluate the sta¬ 
bility of the measured observables against the variation 
of the terms that compose them. For each of the three 
asymmetries extracted, and for each kind of systematic 
effect studied, one of the quantities in the definition of 
the asymmetry was varied two or more times. Then, the 
asymmetries were produced for each variation, for parts 
A and B separately. The two parts were merged as de¬ 
scribed previously. Finally, the systematic uncertainty 
for each bin was computed as 


and obtained using Eqs. (36) and (37), are very small. 
The standard deviation of the differences between the 


&syst 




(41) 
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Source 

BSA 

TSA 

DSA 

7r u background 

0.005 

0.005 

0.009 

Polarization 

0.003 

0.004 

0.008 

Exclusivity cuts 

0.021 

0.019 

0.027 

Transverse correction 

N.A. 

0.006 

0.006 


TABLE IV. Average systematic uncertainties for each source 
of uncertainty and for each asymmetry type. “Polarization” 
stands for Pb, Pt, PbPt for BSA, TSA and DSA, respectively. 


where A t are the asymmetries corresponding to each vari¬ 
ation i out of n variations performed, while Aq is the stan¬ 
dard asymmetry. Each systematic check was performed 
as follows: 


• Systematics on the exclusivity cuts: the “standard” 
analysis was performed with 3cr cut widths (see 
Section V.3), and four others were done, using cut 
widths of 2.25, 2.5, 3.25 and 3.5 a. The bin-by-bin 
systematics were computed using Eq. (41). This 
source is the biggest contributor to the overall sys¬ 
tematic uncertainty (see Table IV), and it encom¬ 
passes also effects due to variations of the dilution 
factor. 


• Systematics on Pb, P t , PbPt'- the beam-spin asym¬ 
metry was computed two more times, taking two 
different values of Pb : Pb + A(Pi) and Pb — A(P&), 
where A (Pb) is the statistical uncertaint y tha t was 
estimated on this quantity (see Section V.6). An 
equivalent treatment was adopted for the TSA and 
the DSA. There is no major kinematic dependence 
for this source of uncertainty, and its contribution 
to the overall systematic uncertainty is very small 
(see Table IVI. 


• Systematics on epn° background subtraction: the 
three asymmetries (BSA, TSA and DSA) were com¬ 
puted each three times, taking three different val¬ 
ues of the Racc factor (defined in Eq. ( |35] )) that was 
used to compute the epn° background. Specifically, 
the background was computed using the “real” 
value Racc increased and decreased by 30%. There 
is some kinematic dependence for this source of un¬ 
certainty, and its contribution to the overall sys¬ 
tematic uncertainty is smaller than that from the 
exclusivity cuts, but larger than those from the po¬ 
larizations or the dilution factors. 


Table IV reports the averages of each kind of system¬ 
atic uncertainty for three three asymmetries. The biggest 
contribution to the systematic uncertainties comes from 
the exclusivity cuts which contributes to ~ 80% of the 
total systematic error. The total systematic uncertainty 
was computed as a quadratic sum of all contributions. 
For all bins and for the three kinds of asymmetries, the 
statistical uncertainty is bigger than the total systematic 
uncertainty. Both kinds of uncertainties are listed, along 
with the values of the asymmetries, in Table |V| 


V.12.1. Radiative corrections 


Afanasev et al. [21] have computed the radiative cor¬ 
rections for the DVCS and BH processes for CLAS kine¬ 
matics. It was found that, given the strict kinematic cuts 
adopted to select the final state, the undetected radiated 
photon can only have small energies. In this case, there¬ 
fore, the main contribution to the radiative correction 
comes from spin-independent soft-photon emission that 
does not affect the polarization observables - while in¬ 
stead it can affect unpolarized cross sections even up to 
the 20% level. The approximation of negligible contribu¬ 
tion from the radiative corrections to the BSA, TSA and 
DSA, compared to the size of the asymmetries, is valid at 
CLAS kinematics at the 0.1% level [21]. Given the statis¬ 
tical uncertainties and the larger size of other systematic 
effects, it was chosen to neglect this contribution. 


VI. MODELS OF GPDS 


In the following sections the experimental asym¬ 
metries are compared to the predictions of four 
GPD models: Vanderhaeghen-Guichon-Guidal (VGG) 
[23] . Goloskokov-Kroll (GK) [25] , Kumericki-Miiller 
(KMM12) ESI, and Goldstein-Gonzalez-Liuti (GLL) [271 . 

Both the VGG and GK models are based on double 
distributions m m to parametrize the (x, £) dependence 
of the GPDs, and on Regge phenomenology for their t 
dependence. The main differences between these two 
models are in the parametrization of the high-t part of 
the electromagnetic form factors and in the fact that 
the parameters of the GK model are tuned using low-xs 
DVMP data from HERA, which are particularly sensitive 
to gluon and sea-quark GPDs. Therefore, the GK model 
is also suited for gluonic GPDs that are not accounted for 
in VGG. However, given the Xb range of the results pre¬ 
sented here, the description of the valence-quark GPDs 
is sufficient. 

KMM12 is a hybrid model designed for global fit¬ 
ting, in which sea-quark GPDs are modeled in a Mellin- 
Barnes integral representation; valence quarks are mod¬ 
eled in terms of these GPDs on the cross-over line £ = 
x ~ xs/(2 — xb)- The parameters of the model were 
fixed using unpolarized-proton DVCS data from CLAS 
and HERMES, as well as the polarized-proton HERMES 
data. The kinematic range of applicability of this model 

o 2 

is defined by the relation 

Finally, the GGL model provides a diquark model 
based parametrization of GPDs that incorporates Regge 
behavior by introducing a spectral function for the spec¬ 
tator diquark’s invariant mass distribution. The param¬ 
eters of the model were obtained by fitting both deep 
inelastic scattering (DIS) structure functions and the 
recently available flavor-separated nucleon form factors 
data [25] , 













17 


VII. RESULTS 

Hereafter, the results for the three asymmetries are 
presented, discussed and compared to the GPD models 
described in Section m The values of the asymmetries 
for each 4-dimensional bin, along with their uncertainties, 
are listed in Table [v| and in Ref. [30] . 

The harmonic structure in <j) of the asymmetries versus 
( Q 2 -xb) and —t was studied by fitting their </> distribu¬ 
tions. For this goal, the fact that the three asymmetries 
have the same denominator has been exploited. In fact, 
by performing a simultaneaous fit of the 3 asymmetries, 
the common denominator can be constrained to be the 
same for the three different observables. Thus BSA, TSA 
and DSA were fitted, respectively, with three functions 
that shared a common denominator (1 + /3 cos (f>). In the 
highest —t bin of the third ( Q 2 -xb ) bin, [3 was set to zero 
due to the limited <fi coverage. 


VII.1. Beam-spin asymmetry 

Figure [l7| shows the beam-spin asymmetry as a func¬ 
tion of 4> for each slice in the Q 2 -x b space (rows) and for 
each bin in —t (columns). Each asymmetry is fitted with 
the function 


c*l u sin </> 
(1 + j3 cos (j>) 


(42) 


and shows a clear sin^-like modulation, with a decreas¬ 
ing amplitude as —t increases, ranging from ~ 0.25 down 
to ~ 0. The dependence in the other kinematic variables 
appears less marked, although a slight decrease in the —t 
slope seems to happen at the highest Q 2 and Xb values. 
This is confirmed by Fig. [l8j which shows the beam-spin 
asymmetry at 90° (i.e. the olu coefficient in Eq. ©) 
as a function of —t for each Q 2 ~Xb bin. The choice of 
the fitting function was motivated by the physics (see 
Eq. @). The data are compared to the predictions of 
the VGG, GK, KMM12 and GGL models. As expected, 
the VGG and GK models don’t show strong differences 
between each other. With respect to the data, they over¬ 
estimate the amplitude of the experimental asymmetries 
especially at low —t and at low Q 2 -Xb ■ At the highest —t 
values, the VGG model gets closer to the data, while the 
GK model is systematically higher. Both models expect 
a steeper —t slope than the data display. This can be 
due to the fact that these models are based on Double- 
Distributions, where the —t dependence is factorized with 
respect to the (x, £) dependence. The data, instead, seem 
to point to more complex correlations between the three 
variables. The GGL model is in good agreement with the 
data at low —t especially for the first and third Q 2 -Xb 
bin, while it diverges away from the data in the high-x^ 
bins. The discrepancy observed for larger xb values is 
an indication that using only DIS and form factor data 
one can only provide a loose constraint on the £ depen¬ 
dence of the model. The best fit to the data is provided 
by the KMM12 model, which however does not cover our 


whole set of kinematics. For consistency, our beam-spin 
asymmetries were also compared to those obtained from 
previous CLAS data [T5] (el-dvcs experiment). For this 
task, the results for the olu coefficient were used, tak¬ 
ing the kinematic bins from the el-dvcs data that were 
closest to our own. The comparison is shown in Fig. [l8j 
where the el-dvcs results are represented by the (green 
online) triangles. The agreement is good, especially con¬ 
sidering the imperfect kinematical overlap. 


VII.2. Target-spin asymmetry 


The results for the target-spin asymmetry m are pre¬ 
sented in Fig. [19] as a function of (f> for each slice in the 
Q 2 -Xb space (rows) and for each bin in —t (columns). As 
for the BSA, it is fitted with the function 


ccul sin cj> 
1 + /3 cos (j) 


(43) 


and shows the typical sin^-like dependence, with ampli¬ 
tudes ranging from 0.1 to 0.3, but its evolution with — t is 
quite different from the BSA, in shape and magnitude. In 
fact, the amplitude of the target-spin asymmetry seems 
rather constant as a function of all kinematic variables, 
—t included, apart from the expected systematic drop 
towards t ~ t rnin . Figure [20] shows the t-dependence 
for each bin in Q 2 -x b of the sin <f> fitting coefficient a\jh 
(Eq. @), which appears fairly constant, unlike what 
was observed for the beam-spin asymmetry. As men¬ 
tioned above, the variable t yields the Fourier conjugate 
of the impact parameter, describing the transverse posi¬ 
tion of the partons in the reference frame where the pro¬ 
ton goes at the speed of light. Therefore, a steep t slope 
is equivalent to a rather flat spatial distribution, and vice 
versa. Equations © and ( fl9] ) point to the proportion¬ 


ality between, respectively, TSA and SmK and BSA and 
SmU. Thus, the t behavior of the TSA compared to 
that of the BSA suggests that the axial charge (linked to 
H) is more concentrated in the center of the proton than 
the electric charge (linked to H). This fact was already 
observed in a paper ;32] devoted to the extraction of the 


CFFs 3mH and An iH from the HERMES data. This 
finding is clearly not predicted by the VGG or GK mod¬ 
els, which instead display a similar drop with t for the 
TSA as what was computed for the BSA. These mod¬ 
els approximately reproduce the low-f magnitude of the 
asymmetry in some kinematics (namely, in Q 2 -Xb bins 
1 and 3), with a slightly better fit of the data for VGG. 
GK predicts an increase of the TSA with xb that is not 
observed in the experiment. By comparing the behavior 
of these two models for the two single-spin asymmetries, 
it can be observed how both reach good agreement with 
the data at high — t for the BSA and a low —t for the 
TSA. These data can therefore provide strong guidance 
to correct the t dependence of the parametrizations of 
both H and H. The GGL model also predicts a drop 
in —t not confirmed by the data, and moreover it over¬ 
estimates the magnitude of the asymmetry by at least a 
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FIG. 17. (Color online) Beam-spin asymmetry for the reaction ep as a function of <j> for the various Q 2 -xb (rows) and 

—t (columns) bins. The point-by-point systematic uncertainties are represented by the shaded bands. The solid black curve is 
the fit with the function in Eq. (421. In the highest —t bin of the third ( Q 2 -xb ) bin, /3 was set to zero due to the limited cj) 
coverage, while no fit is performed on the first —f bin of the highest ( Q 2 -Xb ) bin, where only one data point is present. The 
curves show the predictions of the VGG [23] (red-dashed) and KMM12 [26] (cyan-dotted) models. 
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FIG. 18. (Color online) t dependence for each Q 2 -xb bin 
of the «lu term of the beam-spin asymmetry. The curves 
show the predictions of four GPD models for the BSA at 
<f> = 90°: i) VGG [23) (red dashed), ii) KMM12 |26] (blue 
dotted), iii) GK [26] (blue dash-dotted), and iv) GGL (27] 
(orange dashed-three-dotted). The square black points are 
the results obtained from the present analysis, the triangular 
green data come from the previous CLAS experiment with 
unpolarized proton target |15| . 


factor of 2. The best fit to the data is provided, in the 
bins where it applies, by the KMM12 model. 

Our target-spin asymmetries were compared to those 
obtained from HERMES and to the results of a previous 
exploratory CLAS experiment na. Given the different 
kinematical coverages of the three experiments, it was 
decided to integrate the data over all values of Q 2 -Xb 
of the phase space, and obtain a set of TSAs only as a 
function of (f> and —f, with a finer binning in —t than for 
our four-fold asymmetries. For this comparison, each (f> 
distribution for a given —t bin was fitted with the func¬ 
tion a(j L sin cf> + 7ul sin(2</>), which was the one used for 
both the HERMES and the old CLAS data. The com¬ 
parison of the a{j L coefficients is shown in Fig. [2TJ Our 
data have at least a factor 5 smaller error bars than the 
previously published data, and extend the —t range up 
to 1.7 (GeV/c) 2 . 

A further test was made including in the fit an addi¬ 


tional sin(2</>) term for the TSA. Thus the adopted fitting 
function was 

auL sin (f> + y UL sin(2</>) 

1 + /3 cos </> 

The data from HERMES |TB], in fact, show a non- 
negligible contribution of the sin 2 (f> moment to the TSA. 
It must be pointed out, however, that in the expansion 
in sin <fi moments done by the HERMES Collaboration, 
the contribution of the denominator could mix with the 
sin n(f> terms, while in our analysis it is treated with its 
own parameter. In the Belitsky-Muller-Kirchner formal¬ 
ism, the parts of the DVCS/BH interference term de¬ 
pending on sin 2(f) appear only at twist-three level [9] . 
Therefore, observing a sizeable sin 2 (f> component in the 
target-spin asymmetry would mean having some sensi¬ 
tivity to twist-three CFFs. A recent paper pointed 
to a possible way to access the quarks’ orbital angular 
momentum via the measurement of twist-three GPDs. 
First of all, the stability of Oul was verified by compar¬ 
ing the sin cf> parameter obtained with and without sin 2 (f> 
term in the fitting function. The sin 2(f) term appears 
to be much smaller, at least in the low-f region, than 
the sin <^> term, often compatible with zero, with a slight 
tendency to increase at high —t in some kinematic bins 
towards negative values. The sin </> component is always 
dominant. However, given the limited statistics and the 
relatively small number of (f> bins, the uncertaintes on 
the denominator parameter highly affect the extraction 
of such a potentially small higher-twist sin 2(f> modulation 
in the numerator, so no reliable sin 2(^>-extraction can be 
achieved with the present fitting procedure. 


VII.3. Double-spin asymmetry 


The double-spin asymmetry is plotted in Fig. [22] as a 
function of </> for each bin in —t (columns) and for each 
slice in the Q 2 -Xb space (rows). It is larger in magnitude 
than the single-spin asymmetries presented in the previ¬ 
ous sections (around 0.6), seems rather flat as a function 
of (f>, and presents a slow decrease as a function of — t. 
The data were fitted with the function 


«ll + All cos 4> 
1 + f3 cos (f) 


(45) 


where the denominator parameter /3 is the same as for 
the fits to the two single-spin asymmetries. The two 
sets of fit parameters of the numerator, kll and All, 
are shown, as functions of —t and for each Q 2 -xb bin in 
Figs. [23] and [24] respectively. The constant term dom¬ 
inates the asymmetry, while the cos (f> term of the nu¬ 
merator is compatible with zero for most kinematics. In 


Figs. 22 23 and [24] the three fit parameters of the 
double-spin asymmetries are compared to the four model 
predictions for DVCS+BH and to the calculations for 
BH only (green dot-dashed curve). It seems that Bethe- 
Heitler fully dominates the constant term, and all models 
- except for GGL, which misses both the magnitude and 
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FIG. 19. (Color online) Target-spin asymmetry for the reaction ep —> e'p'j as a function of <j> for the various Q 2 -xb (rows) and 
—t (columns) bins. The point-by-point systematic uncertainties are represented by the shaded bands. The solid black curve is 
the fit with the function in Eq. (431. In the highest —t bin of the third ( Q 2 -xb ) bin, /3 was set to zero due to the limited cj) 
coverage, while no fit is performed on the first —f bin of the highest ( Q 2 -Xb ) bin, where only one data point is present. The 
curves show the predictions of the VGG [23] (red-dashed) and KMM12 [26] (blue-dotted) models. 
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FIG. 20. (Color online) t dependence, for each Q 2 -xb bin, of 
the «ul term of the target-spin asymmetry. The curves show 
the predictions of four GPD models for the TSA at <j> = 90°: 
i) VGG [23] (red dashed), ii) KMM12 [26| (cyan dotted), iii) 
GK [26] (blue dash-dotted), and iv) GGL [27] (orange dashed- 
three-dotted) . 



FIG. 21. (Color online) Comparisons of the t dependences 
of the sine/) term of the ejry target-spin asymmetries for 
the present data, integrated over Q 2 and Xb (black circles), 
the previous CLAS experiment [13] (magenta triangles), and 
HERMES [16] (green squares). 


the t-dependence of this observable - predict this and 
correctly reproduce it. The best match for this term is 
provided by the VGG and GK models, which show size¬ 
able differences only at the highest —t values, where the 
DVCS contribution is expected to start to play a role. 
The models suggest a slight contribution from DVCS in 
the cos (j) term but the statistical precision of the data 
does not allow us to draw conclusions on which predic¬ 
tion provides the better fit. 


VIII. EXTRACTION OF COMPTON FORM 
FACTORS 


In recent years, various groups have developed and 
applied different procedures to extract Compton Form 
Factors from DVCS observables. The approach adopted 
here [34li36] is based on a local-fitting method at each 
given experimental ( Q 2 ,xb , —t) kinematic point. In this 
framework, instead of four complex CFFs defined as in 
Eq. [9j there are eight real CFFs defined as 

F Re (£,t)=$teF(Z,t) (46) 

Fim(€,t) = --$SmF(£,t) = [F(£,£,t) =F F(-£,£,t )], 

7r 

(47) 

where the sign convention is the same as for Eq. (|8|. 
These CFFs are the almost-free parameters - their values 
are allowed to vary within ±5 times the values predicted 
by the VGG model - that are extracted from DVCS ob¬ 
servables using the well-established DVCS+BH theoret¬ 
ical amplitude. The BH amplitude is calculated exactly 
while the DVCS amplitude is taken at the QCD leading 
twist. The expression of these amplitudes can be found, 
for instance, in [25] , 

The three sets of asymmetries (BSA, TSA and DSA) 
for all kinematic bins were processed using this fitting 
procedure to extract the Compton Form Factors. In the 
adopted version of the fitter code, Ej m is set to zero, as 
E is assumed to be purely real - it is parametrized in the 
VGG model by the pion pole (1 /{t — m 2 )). Thus seven 
out of the eight real and imaginary parts of the CFFs are 
left as free parameters in the fit. Figure [25| shows Hj m 
(black full squares) and Hi m (red full circles), which are 
obtained from the fit of the present data, as a function 
of —t for each of our 5 Q 2 -x r bins. These are the two 
CFFs that appear to be better constrained by the present 
results. Given that the size of the error bars reflects 
the sensitivity of the combination of observables to each 
CFF, it is evident that, as expected, our asymmetries are 
mostly sensitive to GmFf. 

The results for _Hj m and F[j m confirm what had been 
previously observed in a qualitative way by direct com¬ 
parison of the f-dependence of our TSAs and BSAs in 
Section |VII.2| the t-slope of SmK is much steeper than 
that of 3m' H, hinting at the fact that the axial charge 
(linked to QnrTf) might be more “concentrated” in the 
center of the nucleon than the electric charge (linked to 
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FIG. 22. (Color online) Double-spin asymmetry for the reaction ep —» e'p'j as a function of 4> for the various Q 2 -xb (rows) and 
—t (columns) bins. The point-by-point systematic uncertainties are represented by the shaded bands. The solid black curve is 
the fit with the function in Eq. (451. In the highest —t bin of the third ( Q 2 -xb ) bin /3, was set to zero due to the limited cj) 
coverage, while no fit is performed on the first —f bin of the highest ( Q 2 -Xb ) bin, where only one data point is present. The 
red-dashed and cyan-dotted curves are predictions of the VGG and KMM12 models, respectively. The pink two-dot-dashed 
curves are the calculations for the Bethe-Heitler process. 
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FIG. 23. (Color online) t dependence for each Q 2 -Xb bin of 
the constant term k ll of the double-spin asymmetry. The 
pink two-dot-dashed curves are the calculations of the DSA 
for the Bethe-Heitler process alone. The curves show the pre¬ 
dictions for the full ejry amplitude of four GPD models: i) 
VGG [22. (red dashed), ii) KMM12 [25] (cyan dotted), iii) 
GK [26] (blue dash-dotted), and iv) GGL [22 (orange dashed- 
three-dotted) . 


FIG. 24. (Color online) t dependence for each Q 2 -xb bin of 
the cos 4> term All of the double-spin asymmetry. The pink 
two-dot-dashed curves are the calculations of the DSA for the 
Bethe-Heitler process alone. The curves show the predictions 
for the full epj amplitude of four GPD models: i) VGG [23] 
(red dashed), ii) KMM12 [25] (cyan dotted), iii) GK [26j (blue 
dash-dotted), and iv) GGL [27] (orange dashed-three-dotted). 


3m?I). This effect seems stronger at the lowest values of 
xb, while both CFFs tend to flatten out as xb increases. 


It is also interesting to compare the results obtained 
for the two equal-a ;b bins ( Q 2 = 1.97 (GeV/c) 2 and 
Q 2 = 2.41 (GeV/c) 2 ): within the limits imposed by the 
size of the error bars and by the Q 2 lever arm (only 0.44 
(GeV/c) 2 ), both sets of CFFs are compatible, at the 1 -a 
level, which supports the validity of the scaling hypoth¬ 
esis. 
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the values of Hj m and Hi m that were ob- 
using the same fitting code with the results 
from Il5l for the beam-spin asymmetry and from [13] for 
the target-spin asymmetry, are also shown. Aside from 
the much larger kinematic coverage for the polarized- 
target observables made accessible by our data, in the 
kinematics where the previous extraction had been at¬ 
tempted, our data improve the precision of SmU. 


IX. CONCLUSIONS 


For the first time four-dimensional single-beam, single¬ 
target, and double (beam-target) spin asymmetries for 
deeply virtual Compton scattering were extracted over 
a large phase space at the same kinematics. This exper¬ 
iment used the CLAS detector in conjunction with the 
IC calorimeter and the Hall-B longitudinally polarized 
14 NH 3 target. 165 4-dinrensional bins in Q 2 , xb , —t 
and </>, covering a wide kinematic range (1 < Q 2 < 5.2 
(GeV/c) 2 , 0.12 < x B < 0.6, 0.08 < -t < 2 (GeV/c) 2 , 
0 < (f> < 360°) were obtained for the three asymme¬ 
tries, with systematic uncertainties largely below the 
statistical uncertainties. The 4> dependence of the 
obtained asymmetries was studied. The dominance 
of the leading-order handbag mechanism is supported 
by the prevalence, especially at low — t, of the sin^> 
term over higher sin ncj) components in both single-spin 
asymmetries. The DSA is found to be mostly dominated 
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-t (GeV 2 ) 


FIG. 25. (Color online) t dependence for each Q 2 -xb bin of 
Him (black squares) and Him (red circles). The full points are 
obtained by fitting the present data (TSA, BSA and DSA). 
The empty points were obtained by fitting the BSA results 
from EH integrated over all values of Q 2 at xb ~ 0.25, and 
the TSAs from [T3j. 
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by the constant term, which contains both BH and 
DVCS/BH interference, although in our kinematics 
the BH contribution is the strongest. These data 
bring important constraints to GPD parametrizations, 
especially for H and H. Using one method among the 
various ones that are currently in development, these 
data allow us to extract the imaginary parts of the H 
and T~L Compton Form Factors and to gain insight, via 
their relative t slopes, about the spatial distribution 
of the electric and axial charges in the proton. The 
extraction of the Compton Form Factors will be further 
improved once the new CLAS results on DVCS cross 
sections will become available [38] . Furthermore, the 
extensive DVCS-devoted program planned at Jefferson 
Lab for the 12-GeV era will extend our knowledge of the 
Generalized Parton Distributions of the proton in terms 
of both kinematical coverage and statistical precision. 

We thank the staff of the Accelerator and the Physics 
Divisions and of the Target Group at Jefferson Lab 
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- 0.0160 

0.63 ± 0.11 ± 0.03 

0.0019 

2.42 

0.260 

0.23 

270 

- 0.23 ± 0 . 06 ± 0.03 

- 0.14 ± 0.08 ± 0.02 

- 0.0118 

0.61 ± 0.09 ± 0.04 

0.0079 

2.48 

0.266 

0.23 

308 

- 0.15 ± 0 . 05 ± 0.02 

- 0.21 ± 0.06 ± 0.01 

- 0.0039 

0.59 ± 0.07 ± 0.02 

0.0115 

2.56 

0.279 

0.22 

336 

- 0.01 ± 0 . 07 ± 0.01 

- 0.09 ± 0.08 ± 0.02 

0.0001 

0.49 ± 0.10 ± 0.03 

0.0127 

2.60 

0.284 

0.41 

31 

- 0.18 ± 0 . 22 ± 0.09 

- 0.01 ± 0.23 ± 0.06 

0.0000 

0.01 ± 0.31 ± 0.11 

- 0.0045 

2.55 

0.275 

0.44 

52 

0.17 ± 0 . 06 ± 0.01 

0.07 ± 0.07 ± 0.02 

0.0055 

0.54 ± 0.08 ± 0.05 

0.0003 

2.45 

0.264 

0.44 

91 

0.20 ± 0 . 07 ± 0.02 

0.04 ± 0.08 ± 0.04 

0.0157 

0.59 ± 0.09 ± 0.03 

0.0022 

2.37 

0.254 

0.44 

122 

0.23 ± 0 . 10 ± 0.03 

0.15 ± 0.12 ± 0.01 

0.0188 

0.69 ± 0.13 ± 0.03 

0.0043 

2.33 

0.250 

0.41 

159 

- 0.15 ± 0 . 13 ± 0.04 

- 0.24 ± 0.16 ± 0.04 

0.0063 

0.59 ± 0.18 ± 0.04 

0.0043 

2.36 

0.253 

0.41 

199 

- 0.02 ± 0 . 12 ± 0.04 

0.00 ± 0.14 ± 0.03 

- 0.0079 

0.59 ± 0.16 ± 0.04 

0.0044 

2.36 

0.254 

0.43 

235 

- 0.28 ± 0 . 10 ± 0.03 

- 0.05 ± 0.12 ± 0.03 

- 0.0184 

0.71 ± 0.13 ± 0.04 

0.0045 

2.43 

0.261 

0.45 

270 

- 0.26 ± 0 . 07 ± 0.02 

- 0.12 ± 0.09 ± 0.03 

- 0.0161 

0.60 ± 0.10 ± 0.02 

0.0023 

2.54 

0.274 

0.45 

307 

- 0.24 ± 0 . 07 ± 0.03 

- 0.22 ± 0.08 ± 0.03 

- 0.0065 

0.43 ± 0.09 ± 0.05 

- 0.0004 

2.61 

0.284 

0.42 

330 

- 0.15 ± 0 . 19 ± 0.07 

- 0.03 ± 0.18 ± 0.03 

0.0004 

0.59 ± 0.27 ± 0.14 

0.0000 

2.43 

0.261 

1.22 

20 

- 0.01 ± 0 . 08 ± 0.02 

0.08 ± 0.09 ± 0.01 

0.0007 

0.30 ± 0.11 ± 0.02 

- 0.0068 

2.45 

0.265 

1.06 

52 

0.14 ± 0 . 09 ± 0.03 

- 0.00 ± 0.11 ± 0.03 

0.0034 

0.60 ± 0.12 ± 0.02 

- 0.0035 

2.37 

0.255 

0.96 

88 

0.18 ± 0 . 20 ± 0.12 

0.21 ± 0.22 ± 0.05 

0.0085 

0.65 ± 0.26 ± 0.07 

0.0004 

2.43 

0.261 

0.96 

273 

- 0.31 ± 0 . 17 ± 0.05 

0.09 ± 0.20 ± 0.08 

- 0.0069 

0.49 ± 0.23 ± 0.06 

0.0006 

2.44 

0.262 

1.08 

308 

- 0.16 ± 0 . 10 ± 0.04 

- 0.04 ± 0.11 ± 0.03 

- 0.0037 

0.45 ± 0.13 ± 0.05 

0.0045 

2.46 

0.265 

1.23 

340 

- 0.04 ± 0 . 07 ± 0.02 

- 0.09 ± 0.09 ± 0.01 

- 0.0007 

0.44 ± 0.10 ± 0.03 

0.0059 

2.55 

0.319 

0.16 

19 

0.13 ± 0 . 05 ± 0.02 

0.07 ± 0.07 ± 0.01 

0.0020 

0.73 ± 0.07 ± 0.06 

0.0546 

2.50 

0.319 

0.16 

50 

0.25 ± 0 . 07 ± 0.03 

0.10 ± 0.09 ± 0.01 

0.0045 

0.60 ± 0.10 ± 0.06 

0.0513 

2.53 

0.320 

0.16 

88 

0.37 ± 0 . 14 ± 0.08 

0.16 ± 0.17 ± 0.03 

0.0066 

0.74 ± 0.19 ± 0.05 

0.0464 

2.55 

0.317 

0.16 

240 

- 0.13 ± 0 . 22 ± 0.01 

- 0.19 ± 0.26 ± 0.05 

- 0.0064 

0.55 ± 0.29 ± 0.03 

0.0345 

2.59 

0.321 

0.16 

273 

- 0.12 ± 0 . 15 ± 0.05 

- 0.24 ± 0.17 ± 0.07 

- 0.0078 

0.93 ± 0.19 ± 0.09 

0.0479 

2.53 

0.319 

0.16 

309 

- 0.21 ± 0 . 07 ± 0.02 

- 0.11 ± 0.08 ± 0.01 

- 0.0045 

0.77 ± 0.09 ± 0.04 

0.0528 

2.55 

0.318 

0.16 

342 

- 0.12 ± 0 . 06 ± 0.02 

0.02 ± 0.07 ± 0.01 

- 0.0005 

0.70 ± 0.08 ± 0.05 

0.0534 

2.72 

0.338 

0.24 

19 

0.08 ± 0 . 03 ± 0.00 

0.10 ± 0.03 ± 0.01 

0.0000 

0.72 ± 0.04 ± 0.03 

0.0398 

2.68 

0.338 

0.24 

51 

0.20 ± 0 . 03 ± 0.01 

0.11 ± 0.04 ± 0.01 

0.0050 

0.69 ± 0.05 ± 0.03 

0.0387 

2.67 

0.337 

0.24 

89 

0.16 ± 0 . 05 ± 0.02 

0.15 ± 0.07 ± 0.02 

0.0223 

0.69 ± 0.07 ± 0.03 

0.0293 

2.43 

0.337 

0.25 

122 

0.13 ± 0 . 07 ± 0.01 

0.19 ± 0.08 ± 0.01 

0.0373 

0.72 ± 0.09 ± 0.04 

0.0164 

2.47 

0.326 

0.25 

161 

0.14 ± 0 . 12 ± 0.01 

0.25 ± 0.14 ± 0.03 

0.0214 

0.59 ± 0.15 ± 0.07 

- 0.0042 

2.55 

0.322 

0.25 

201 

- 0.07 ± 0 . 09 ± 0.03 

- 0.25 ± 0.11 ± 0.06 

- 0.0179 

0.75 ± 0.12 ± 0.02 

- 0.0020 

2.55 

0.338 

0.24 

238 

- 0.20 ± 0 . 07 ± 0.02 

- 0.18 ± 0.09 ± 0.01 

- 0.0346 

0.71 ± 0.09 ± 0.03 

0.0142 

2.73 

0.337 

0.24 

271 

- 0.26 ± 0 . 06 ± 0.03 

- 0.28 ± 0.07 ± 0.02 

- 0.0233 

0.68 ± 0.07 ± 0.03 

0.0278 

2.69 

0.336 

0.24 

309 

- 0.26 ± 0 . 03 ± 0.02 

- 0.16 ± 0.04 ± 0.01 

- 0.0058 

0.73 ± 0.05 ± 0.03 

0.0385 

2.68 

0.338 

0.24 

341 

- 0.04 ± 0 . 03 ± 0.00 

- 0.08 ± 0.03 ± 0.01 

0.0002 

0.59 ± 0.04 ± 0.03 

0.0395 
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o z 

(GeV/c) 1 2 3 4 5 6 

XB 

— t 

(GeV/c) 2 

0 

(deg) 

BSA ± stat. ± syst. 

TSA ± stat. ± syst. 

CAut 

DSA ± stat. ± syst. 

Ca lt 

2.81 

0.356 

0.44 

20 

0.09 ± 0.03± 0.01 

0.04 ± 0.03 ± 0.01 

-0.0055 

0.59 ± 0.04 ± 0.02 

0.0123 

2.74 

0.346 

0.46 

51 

0.20 ± 0.03± 0.01 

0.13 ± 0.03 ± 0.01 

0.0049 

0.66 ± 0.04 ± 0.02 

0.0143 

2.68 

0.346 

0.45 

90 

0.15 ± 0.05± 0.02 

0.18 ± 0.06 ± 0.02 

0.0358 

0.68 ± 0.06 ± 0.02 

0.0135 

2.24 

0.349 

0.47 

126 

0.07 ± 0.05± 0.01 

0.19 ± 0.06 ± 0.03 

0.0636 

0.53 ± 0.06 ± 0.02 

0.0097 

2.25 

0.339 

0.49 

162 

0.06 ± 0.06± 0.01 

0.25 ± 0.07 ± 0.03 

0.0325 

0.68 ± 0.08 ± 0.02 

0.0061 

2.23 

0.339 

0.49 

202 

-0.13 ± 0.06± 0.02 

-0.15 ± 0.07 ± 0.02 

-0.0296 

0.54 ± 0.08 ± 0.04 

0.0022 

2.43 

0.347 

0.47 

232 

-0.06 ± 0.06± 0.01 

-0.21 ± 0.08 ± 0.03 

-0.0567 

0.38 ± 0.09 ± 0.02 

0.0034 

2.72 

0.344 

0.47 

271 

-0.22 ± 0.05± 0.02 

-0.26 ± 0.06 ± 0.02 

-0.0362 

0.62 ± 0.06 ± 0.04 

0.0118 

2.77 

0.347 

0.45 

309 

-0.21 ± 0.03± 0.01 

-0.18 ± 0.04 ± 0.01 

-0.0058 

0.64 ± 0.04 ± 0.02 

0.0137 

2.77 

0.355 

0.45 

340 

-0.15 ± 0.02± 0.01 

-0.07 ± 0.03 ± 0.01 

0.0049 

0.59 ± 0.03 ± 0.02 

0.0129 

2.59 

0.352 

1.30 

20 

0.02 ± 0.03± 0.00 

0.04 ± 0.04 ± 0.01 

-0.0041 

0.36 ± 0.04 ± 0.03 

-0.0051 

2.59 

0.348 

1.09 

52 

0.13 ± 0.04± 0.01 

0.15 ± 0.05 ± 0.02 

0.0088 

0.39 ± 0.06 ± 0.02 

-0.0035 

2.28 

0.350 

1.28 

87 

0.12 ± 0.06± 0.02 

0.09 ± 0.08 ± 0.06 

0.0326 

0.36 ± 0.08 ± 0.03 

-0.0008 

2.27 

0.353 

1.09 

129 

0.08 ± 0.06± 0.03 

0.19 ± 0.08 ± 0.05 

0.0375 

0.30 ± 0.08 ± 0.04 

-0.0025 

2.33 

0.345 

1.06 

161 

-0.00 ± 0.06± 0.02 

0.21 ± 0.08 ± 0.03 

0.0188 

0.52 ± 0.08 ± 0.02 

0.0015 

2.29 

0.347 

1.06 

201 

-0.00 ± 0.06± 0.02 

-0.19 ± 0.08 ± 0.03 

-0.0184 

0.52 ± 0.08 ± 0.02 

0.0020 

2.33 

0.352 

1.23 

236 

-0.01 ± 0.06± 0.02 

-0.25 ± 0.09 ± 0.06 

-0.0384 

0.40 ± 0.09 ± 0.02 

0.0003 

2.25 

0.347 

1.23 

269 

-0.04 ± 0.05± 0.01 

-0.19 ± 0.06 ± 0.04 

-0.0356 

0.30 ± 0.06 ± 0.02 

-0.0027 

2.61 

0.349 

1.09 

308 

-0.11 ± 0.05± 0.01 

-0.14 ± 0.06 ± 0.02 

-0.0086 

0.48 ± 0.06 ± 0.03 

-0.0016 

2.57 

0.352 

1.28 

340 

-0.02 ± 0.03± 0.00 

-0.07 ± 0.03 ± 0.01 

0.0036 

0.43 ± 0.04 ± 0.02 

-0.0034 

3.21 

0.414 

0.27 

19 

-0.24 ± 0.18± 0.02 

0.00 ± 0.00 ± 0.00 

0.0000 

0.90 ± 0.21 ± 0.02 

0.0046 

3.16 

0.414 

0.27 

51 

0.22 ± 0.21± 0.09 

-0.04 ± 0.26 ± 0.08 

-0.0007 

0.30 ± 0.29 ± 0.09 

0.0070 

3.36 

0.443 

0.51 

17 

0.07 ± 0.04± 0.01 

-0.02 ± 0.04 ± 0.01 

-0.0034 

0.82 ± 0.05 ± 0.03 

0.0398 

3.41 

0.444 

0.49 

52 

0.21 ± 0.05± 0.02 

0.13 ± 0.06 ± 0.02 

0.0066 

0.77 ± 0.06 ± 0.04 

0.0395 

3.25 

0.448 

0.49 

89 

0.24 ± 0.07± 0.03 

0.27 ± 0.09 ± 0.04 

0.0035 

0.77 ± 0.10 ± 0.05 

0.0397 

2.84 

0.430 

0.53 

124 

0.12 ± 0.09± 0.05 

0.27 ± 0.11 ± 0.03 

0.0548 

0.64 ± 0.12 ± 0.07 

0.0284 

2.93 

0.420 

0.57 

155 

-0.08 ± 0.24± 0.08 

0.07 ± 0.31 ± 0.08 

0.0227 

0.44 ± 0.30 ± 0.08 

0.0040 

2.94 

0.421 

0.58 

207 

-0.44 ± 0.23± 0.06 

-0.25 ± 0.28 ± 0.07 

-0.0273 

0.50 ± 0.31 ± 0.07 

0.0082 

2.95 

0.435 

0.50 

238 

-0.04 ± 0.11± 0.05 

-0.18 ± 0.13 ± 0.05 

-0.0497 

0.75 ± 0.14 ± 0.07 

0.0309 

3.33 

0.440 

0.48 

271 

-0.08 ± 0.08± 0.02 

-0.24 ± 0.10 ± 0.01 

-0.0322 

0.82 ± 0.11 ± 0.04 

0.0381 

3.46 

0.445 

0.50 

310 

-0.28 ± 0.05± 0.02 

-0.23 ± 0.06 ± 0.01 

-0.0879 

0.68 ± 0.06 ± 0.04 

0.0364 

3.35 

0.443 

0.52 

342 

-0.08 ± 0.03± 0.01 

-0.16 ± 0.04 ± 0.01 

-0.0010 

0.77 ± 0.04 ± 0.03 

0.0392 

3.37 

0.465 

1.12 

18 

0.06 ± 0.03± 0.01 

0.08 ± 0.04 ± 0.01 

-0.0064 

0.51 ± 0.04 ± 0.02 

0.0015 

3.46 

0.465 

1.13 

52 

0.21 ± 0.05± 0.02 

0.16 ± 0.07 ± 0.03 

0.0095 

0.60 ± 0.07 ± 0.03 

0.0071 

3.18 

0.459 

1.27 

89 

-0.00 ± 0.07± 0.02 

0.15 ± 0.10 ± 0.05 

0.0507 

0.55 ± 0.10 ± 0.04 

0.0078 

3.13 

0.456 

1.15 

126 

0.11 ± 0.06± 0.01 

0.04 ± 0.08 ± 0.05 

0.0558 

0.55 ± 0.08 ± 0.02 

0.0007 

3.30 

0.434 

1.18 

159 

0.07 ± 0.11± 0.09 

0.12 ± 0.15 ± 0.04 

0.0241 

0.43 ± 0.15 ± 0.07 

-0.0116 

3.11 

0.430 

1.09 

203 

-0.10 ± 0.13± 0.03 

0.00 ± 0.16 ± 0.06 

-0.0219 

0.78 ± 0.16 ± 0.06 

-0.0021 

3.16 

0.455 

1.25 

236 

-0.06 ± 0.07± 0.01 

-0.30 ± 0.10 ± 0.06 

-0.0644 

0.46 ± 0.09 ± 0.03 

-0.0026 

3.12 

0.454 

1.35 

269 

-0.02 ± 0.06± 0.02 

-0.21 ± 0.09 ± 0.05 

-0.0541 

0.56 ± 0.08 ± 0.03 

0.0083 

3.49 

0.467 

1.10 

310 

-0.31 ± 0.06± 0.04 

-0.25 ± 0.07 ± 0.03 

-0.0119 

0.77 ± 0.08 ± 0.05 

0.0115 

3.36 

0.464 

1.15 

342 

-0.07 ± 0.03± 0.01 

-0.08 ± 0.04 ± 0.01 

0.0063 

0.51 ± 0.04 ± 0.02 

0.0014 


TABLE V: Values of the three asymmetries with their statistical and 
systematic uncertainties for each 4-dimensional bin. The values of the 
bin-by-bin transverse corrections for the TSA and DSA are also shown. 
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